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Abstract. We present easy to verify conditions implying stability estimates 
for operator matrix splittings which ensure convergence of the associated Trot- 
ter, Strang and weighted product formulas. The results are applied to inho- 
mogeneous abstract Cauchy problems and to boundary feedback systems. 



1. Introduction 

Many systems in physics, biology or engineering can be described by an abstract 
Cauchy problem of the form 



(ACP) 



-ftU{t) = Au{t) for t > 0, 
u{0) = Uq 



on a product £ = E x F oi two Banach spaces E and F, see Batkai and Piazzera 
g], Engel and Nagel Chapter VI], Casarino et al. [6j, or Tretter [K]. By Engel, 
Nagel [HI Section II. 6] the problem (|ACP|) is well-posed if and only if the system 
operator A is the generator of a strongly continuous semigroup (7"(i))t>o on £. 
Moreover, in this case the unique (mild) solution u{-) of (lACPp is given by 

However, in general it is not possible to calculate the entries of 

in terms of A in order to obtain an explicit representation of the solution u{-). But 
as we will see below this can be achieved in case A has some special structure, e.g., 
if A is of triangular form. The idea at this point is to split A into (a sum of) simpler 
pieces, for which it is possible to calculate the associated semigroup and then to use 
some kind of product formula to reassemble T{t) from these pieces. This approach 
is made more precise in the following result. 

Theorem 1.1. For z = 1,2 let Ai be the generator of the strongly continuous 
semigroup {Ti{t))t>o on the Banach space £. Suppose that A := Ai + A2 is the 
generator of the strongly continuous semigroup (T(t))t>o- Then the following as- 
sertions are true, 
(i) If there exists M > 1 and a; G M such that 

(1.1) ||(r2(^)ri(^))"|| <Me"* for allt>0 andn&N, 
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then for all x G £ 
(1.2) Tit)x^ lim (r2(^)ri(^))"^, and 

(1-3) rw^ = (ri(^)r2(^)ri(^))" ^. 

(a) If there exists M > I and a; G M such that 



(1.4) 



^{T^{i)T,{r) + T,(im{i)y 



< A/c"* for allt>0 and n e N, 



2" 

t/ien /or all x ^ £ 

(1-5) = „ii,m (ri(^)r2(i) + r2(i)ri(i))" ^. 

For the proofs we refer to Engel and Nagel PI Corollary III. 5. 8], Csomos and Nickel 
[HI Section 2], and Batkai, Csomos and Nickel Section 4]. 

Product formulas like (|1.2I) . (|1.3p and (|1.5I) have been applied to approximate 
the solution of a variety of complicated differential equations and are referred to 
as "operator splitting" in numerical analysis, see for example the monographs by 
Far ago and Havasi [TUl, Holden et al. [T^ or Hundsdorfer and Verwer [T3]. The 
procedure described in Equation (|1.2p is called the Trotter product formula, or 
sequential splitting. Equation (jl.3|) is called the Strang splitting, and Equation 
(|1.5|) is called the (symmetrically) weighted splitting or additive operator splitting. 
These and many other different procedures have been introduced to increase the 
order of convergence. In the finite dimensional setting, sequential splitting is of first 
order, while the other two are of second order. There are many more higher order 
methods in the literature, see Hairer, Lubich and Wanner [TTl Section III. 5. 4], but 
we concentrate here on these three main cases since they are the most frequently 
used ones in applications. 

Various generalizations of this procedure are possible but will not be considered 
in this paper. For non-autonomous versions of these product formulas we refer to 
Batkai et al. [1]. For the combined effect of spatial approximation and operator 
splitting see Batkai, Csomos and Nickel [3], and for the combination of rational 
approximations, operator splitting and spatial approximation see Batkai et al. [5]. 

The crucial hypothesis to achieve convergence of these splitting procedures are 
stability conditions like (|l.ip or (|1.4p . In case the semigroups involved are not 
quasi-contractive, it is in general very difficult to verify these conditions by explicit 
computations. 

The aim of this paper is to address this problem for a special class of triangular 
matrix operator semigroups, which occur quite frequently in applications. To this 
end, in Section [5] we investigate the stability of the Trotter, Strang and weighted 
product formulas for triangular operator matrices. To do that we first characterize 
generators of triangular operator matrix semigroups. Then we analyze the condi- 
tions (|l.ip and (|1.4p in the triangular case and give an abstract sufficient condition 
ensuring them. Finally, we show how extrapolated Favard classes can be used 
to obtain the desired estimates. In Section |3] we consider two classes of applica- 
tions: Inhomogeneous abstract Cauchy problems and abstract boundary feedback 
systems. 

In what follows we use the term "semigroup" to indicate a strongly continuous 
one-parameter semigroup of bounded linear operators, our main reference on this 
topic is Engel, Nagel [9]. 
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2. Splitting for Operator Matrices 

In this section we first characterize generators of triangular matrix semigroups. 
Then we present conditions implying stability for products of triangular operator 
matrix semigroups. Finally, we show how our main assumption on the growth of 
the off-diagonal elements of the matrix semigroup can be verified by the use of 
Favard classes. 

2.1. Characterization of Triangular Matrix Semigroups. As mentioned al- 
ready in the introduction, in general it is not possible to give an explicit matrix 
representation of a semigroup {T{t))t>o on a product space in terms of the en- 
tries of the associated generator. However, things get much simpler if we restrict 
our attention to matrices of triangular form. In order to characterize this class of 
operators we associate to an operator 

A : D{A) C£ 

defined on the product space £ — E x F the operator A : D{A) Q E ^ E defined 

by 

Ax := 7ri(.A(^)) for x £ D{A) := [z E : (^) e D{A)] , 

where tt^ denotes the projection on the i**^ coordinate. Moreover, we denote by 
s{A) the spectral bound of A. With these notations the following result holds. 

Proposition 2.1. Let A generate the semigroup T — {T{t))t>o on £ — E x F . 
Then T has upper triangular form, i.e. 

(2.1) m - (^^^^ for all t > 0, 

if and only if the following two conditions are satisfied. 

(i) For all X e E with (q) G D{A) we have tt2{A{'^)) = 0. 
(ii) There exists X G piA) satisfying Re A > s(A). 

Moreover, in this case (T{t))t>o is a semigroup with generator A. 

Proof. Note first that, if T{t) has upper triangular form (|2.ip . then the entries T{t) 
and S{t) form semigroups. Denote their generators by A and B, respectively. By 
taking the Laplace transform of t M> T{t) we obtain for A large that 

«(A,.4)^(«<y) 

i.e., R{X,A) has upper triangular form for A large. Conversely, if R{X,A) has 
triangular form for sufficiently large A, the Post-Widder inversion formula (see 
Engel and Nagel O Corollary III. 5. 5]) implies that T{t) is upper triangular. Hence 
T{t) is of upper triangular form for all t > if and only if R{X, A) is of upper 
triangular form for all A sufficiently large. This is further equivalent to the fact 
that for some A G C satisfying Re A > s{A) the resolvent has upper triangular 
form. To see this we note that for |A — ^| < ||i?(A, ^)||^^ we have /i G p{A) and 

i?(Ai,^) = 5](A-Ai)"i?(A,^r+'- 

Here the right-hand side yields matrices of upper triangular form and by holomor- 
phy of the resolvent map we conclude that A) is of upper triangular form in 
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the whole connected component of p{A) which is unbounded to the right. After 
these preparations we turn to the proof. 

Suppose that T{t) has upper triangular form for alH > and take some (q) € D{A). 

Then wc obtain 7r2(T(t)(Q) — (q)) =0 and (i) follows by the definition of the 
generator of a semigroup. To show (ii) wc fix A e piA) sufficiently large such that 



R{X,A) 



Ri R2 
R4 



We prove that i?i is the inverse of A — ^, i.e., A G p{A) and i?i — R{X,A) (which 
also implies that A is the generator of {T{t))t>o)- Indeed, for an arbitrary x £ E, 
we have 

R{X,A)0 = {''^^)eD{A), 
and hence by definition Rix £ D{A), i.e., rgi?i C D{A). Moreover, 

(A-^)i?ix = 7ri((A-^)(«^")) =7ri((A-^)i?(A,^)(2)) =a;, 

i.e., i?i is the right-inverse of A — ^. We show that it is also a left-inverse. For 
X e D{A) we have 

i?i(A - A)x = R,n,[{X - = TTi (^R{X,A)(^'''^^^ 'MDY 

which, by validity of (i), further equals to 

= n,[R{X,A){X-A){l))=n,{l)=x. 

Summing up, A £ p{A), hence (ii) is true. Moreover, this implies that A = A. 

Suppose now that (i) and (ii) are satisfied, and fix a A € p{A) n p(A). We have to 
prove that 

i?(A, A) = takes the form ^ 

i.e., i?3 = or, equivalently, 'k-2{R{X,A){^) = for all x G E. Take x G E and 
consider the vector 



R{X,A){1) 



R{X,A)x\ _ /RixX fR{X,A)x 
) ^ [r^xJ "I 



which belongs to ker(A — .4) = {0}. Indeed, we have 
and by (i) 

Hence i?(A, A)x = R\x and iJsa; = 0, and the proof is completed. □ 
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2.2. Stability Conditions for Matrix Products. We recall that the underlying 
idea of our approach is to split a given operator matrix A generating a semigroup 
(7~(t))t>o on a product space £ into a sum A = A\ + Ai of simpler, i.e. triangular, 
matrices Ai^ i = 1,2, and then compute T{t) using some (e.g. the Trotter) product 
formula. Here the crucial hypothesis for convergence is a stability condition on the 
products of the triangular semigroups (7i(t))t>o, see and ()1.4p in Theorem ll.il 
In this section we will consider three types of such splittings and deduce con- 
ditions ensuring that the related stability conditions are satisfied. We start by 
considering two operator matrix semigroups of upper triangular form and ask for 
conditions ensuring that the associated stability condition for the product is satis- 
fied. Let us investigate first the stability condition (|l.ip for the sequential splitting 
(|1.2p and the Strang splitting (|1.3|) . We remark here that the Strang splitting is 
precisely then stable, when the sequential splitting is. Furthermore, the stability 
assumption as in (|l.ip is equivalent to 

||(7^(^)7'i(^))''|| < Me""*^ for a\\ t > and n, fc e N. 

This is trivially true for all splittings considered in this paper and will be used 
without further reference (replace t hy ^ and interchange the roles of n and k). 
The equivalence of the estimates above is even true for more general finite difference 
schemes, the special splitting structure plays no role here. 

Theorem 2.2. Suppose that for i = 1,2 the matrix Ai generates on £ — E x F 
the semigroup (7i(i))t>o of upper triangular form 



T,{t) R,{t) 
S,{t) 



If there exist M' > 1, K > and a;' e M such that for all i = 1,2, t > and n eN 
(^) ||(T2(l)Ti(A))"|| < M'e"'* and ||(52(^)5i(^))"|i < M'e"'S 
(li) \\Ri{t)\\ < ift -e"'*, 

then there are M > I and w G M such that 

(2.2) ||(75(^)Ti(^))"|| < Afe"* for allt>0 andn eE. 

Proof. Since strongly continuous semigroups are exponentially bounded, we can 
choose M' > 1 and cj' G M without loss of generality so that 

mm<M'e^'', and ||5,(t)|| < M'e"'* 

are satisfied for i = 1,2. 

For h > calculate the product 



T2ih)Ti{h) 



T2{h) R2{h)\ fT,{h) Ri{h) 
S2{h)) \ Si{h) 

T2ih)Ti{h) T2{h)Ri{h) + R2{h)Si{h) 
S2{h)Si{h) 



and set R(h) := T2{h)Ri{h) + R2{h)Si{h). This implies 

iT2{h)T,{h)f = i T2{h)T,{h)R{h) + R{h)S2{h)S^{h)\ ^ 
^ ^ I {S2{h)Si{h)f j ' 
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and by induction one can show that 



(2.3) 

where 

(2.4) 



{S2(Ji)S^{K)) 



k 1 



fe-1 



In order to prove (12. 2p . we only have to show the exponential estimate for (★), the 
other entries of the product fulfill such estimates by assumption. Since l|i?(/i)|| < 
2M'Khe'^'^''\ this implies 



fc-i 



{*)\\<Y.\\iT2{hm{h)y -wRm- {s,{h)s,{h)) 



< 2M'3i\:/i^e"'^''e2"'''e"'(''-^"^^'' = 2M'^Khke'^''-''+'^^'' . 
j=o 

If we set h = ^ and k = n we get for M := 2M'^K and cj := w' + |a;'| + 1 

11(^)11 < Mte("'+I"^'l)* < Me"*. 

This completes the proof. 



□ 



In the same spirit and using analogous calculations, we can investigate the stability 
condition (|1.4p for the weighted splitting (|1.5p . 

Theorem 2.3. Suppose that for i — 1,2 the matrix At generates on £ = E x F 
the semigroup {Ti{t))t>o of upper triangular form 



T,{t) R,{t) 
S,{t) 



If there exist M' >l, K > and w' £ M such that for all i ^ 1,2, t >Q and n e M 
(^) \\UTl{i)T2{^)+T2{i)T^{^,)T\\<M 



and 



(ii) \\R,{t)\\ <Kt-e^'\ 
then there are M > 1 and w G R such that 



(2.5) 



< Me'^ 



for allt>{) and n £ N. 



Proof. Again, since strongly continuous semigroups are exponentially bounded, we 
can choose M' > 1 and cj' G R without loss of generality so that 

||r,(t)|| < M'e'^'\ and ||5,(i)|| < M'e"'* 

are satisfied for i = 1,2. 
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Using the computations of the proof of Theorem 12.21 we obtain for h > that 



Ti{h)T2ih) + T2ih)Ti{h) 

'Ti{h)T2{h)+T2{h)Ti{h) M 

Si{h)S2ih) + S2{h)Siih) 

where = Ti{h)R2{h) + Ri{h)S2{h) + T2{h)Ri{K) + R2{h)Si{h). 

Let R'{h) := Ti{h)R2{h) + Ri{h)S2{h) + T2{h)Ri{h) + R2{h)Si{h). Then by 
induction one can verify the identity 



(2.6) {Ti{h)T2{h) + T2{h)Ti{h)y 



{Ti{h)T2{h)+T2{h)Ti{h)y 



(**)/c 



{Si{h)S2{h) + S2(h)Si(h)y 



k > 



where 



(2.7) 



Mfc = ^(Ti(/i)r2(M + r2(/i)ri(/i))'i?'(/i)- 

• (^i(/l)^2(/l) + 5'2(/l)5l(/l)) 



In order to prove (|2.5p . we only have to show the exponential estimate for {**)k, the 
other entries of the product fulfill such estimates by assumption. Since ||i?'(ft.)|| < 
AM'Khe^"''\ this implies 



fc-i 



i**h\\ < J2 \\{Tiih)T2ih) + T2ih)T^ih)y ■ \\R'ih)\\ 

3=0 

■\{S^{h)S2{h) + S2{h)Si{h)f'^'' 



k-1 



= 2M'^Khke'^'^^+^^^2^. 



If we set /i = - and k ^ n wc set for AI := 2M' K and lo :— uj' 

n ° 

II (★★)„!! < Mie("'+I"^'l)*2" < Me"'*2". 
Combining these estimates, the desired statement (|2.5p follows. 



□ 



Summing up, Theorems 12.21 and 12.31 show that the stability condition in (i) for the 
diagonal entries combined with the growth estimate in [ii) imply stability for the 
matrix products. In the next subsection we will come back to condition {ii). 

But first we consider the following stability result for the Trotter, Strang and 
weighted splitting, which does not make use of a special matrix structure. However, 
in Subsection 13.21 we will apply them in the context of matrix decompositions. 



Proposition 2.4. Let A generate a semigroup {T{t))t>o on the Banach space £ 



and denote by {S(t))t>o the semigroup generated by C £ i-c, S{t) = e 
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Then there exist constants M >\ and w G M such that for all t > and n e N 

||(5(A)r(^))"||<Me-* and 



^('5(^)r(^) + r(^)5(l))" 



Proof. By Engel and Nagel [9l Lemma II. 3. 10], there exists an equivalent norm ||| • ||| 
on £ such that (7~(i))t>o is quasi-dissipative for ||| • |||, i.e., satisfies an estimate 

|||T(t)||| < c"'* forallt>0 

and some u' G M. Moreover, 

< e"!'^"!* foralH>0, 

where |||C||| denotes the operator norm of C G JC{£) induced by ||| • |||. Since || • || — ||| ■ ||| 
there exist m', M' > such that m'||A'|| < fxl < M'\\x\\ for all x & £ and hence 

\\imni)r4<ikr-mmi)r4\ 

< ^ . ellicill* . e"'* . ||U||| 

— m III III 

< m; (iiiciiiw)t . ii^ll 

for all t > and x G 8. This implies the first estimate for M := ^ and uj := 

— ^ 'in' 

uj' + |||C|||. The second estimate follows similarly from 

II inmi) + 3im^)r4 < ■ (2 • e«i'^i«- ■ e-'^)" . im 

< 2" • • e(lll'^lll+"')* ■ WxW 
for the same constants M and ut as above. □ 
The previous result applies in particular to the splitting 

Ao + C where C ^ (^^ q 

for some C £ C{E,F), if we assume that Ao generates a matrix semigroup on 
£ = E X F. In this case the semigroup {S{t))t>o generated by C is given by 



(2.8) Sit) = 



Id 
tC Id 



2.3. Estimates for Triangular Matrix Semigroups. As we saw in the previous 
subsection, cf. condition (iii) in Theorems l2.2l and l2.3[ in order to obtain the desired 
stability estimates (12. 2p and (|2.5p we need estimates of the type |ji?i(t)|| < Kt, 
i = 1,2, for the upper right entries Ri{t) of Ti{t). In this section we will use an 
approach based on the concept of Favard classes to achieve this goal. 

Recall from the proof of Proposition l2.1l that given a matrix semigroup (T(i))t>o 
of a triangular form (|2.1I) the diagonal entries {T{t))t>o and {S{t))t>o are semi- 
groups on E and F, respectively. If A and B denote their generators, we define the 
diagonal matrix 

^=(0 b)' DiV) = DiA) X DiB) 
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which generates the diagonal semigroup {S{t))t>o given by 

Sit) 

Moreover, we denote by (5-i(i))t>o the extrapolated semigroup 

S-i{t) 



(Tit) 
I S{t) 



T-i{t) 
S-iit) 



with generator 



^?-i=(^o' /J, D{V^,) = DiA^^) X D{B.^), 

see Engel and Nagel [9l Section II. 5]. For the convenience of the reader we collect 
here some facts concerning Favard classes of semigroup generators. A much more 
detailed account can be found in Engel and Nagel P| Section II. 5. b]. 

Definition 2.5. Let {T{t))t>a be a strongly continuous semigroup on a Banach 
space E with generator A. Then we define its Favard class (or space) as 

Favi(A) -.^ {x €E : sup t~^ ■ \\T{t)x - x\\ < ooj C -B, 

^ te(o,i] 

which becomes a Banach space with respect to the norm 

||a;||Favi(A) := Ikll + sup t''^ ■ \\T{t)x - x\\. 

te{o,i] 

We note that for reflexive Banach spaces E one always has Favi(A) = D{A) (see 
Engel and Nagel ^ Corollary II. 5. 21]), hence Favard spaces are interesting only in 
nonreflexive spaces. 

One can define the Favard space Favo(A) = Favi(A_i) for the extrapolated semi- 
group (T-i(t))t>o with generator A_i in a similar manner. Using these notations 
we have the following result. 

Proposition 2.6. Let (7~(t))t>o be a triangular semigroup of the form (|2.ip on 
the product space £ = E x F with generator A. Then the following assertions are 
equivalent. 

(a) There exists K >0 such that \\R{t)\\ < Kt for all t e [0, 1]. 

(a') There exists K >0, lo eR such that \\R{t)\\ < Kt ■ e"* for all t > 0. 

(b) There exists P £ C{F,Y&yo{A)) such that A = (2?-i + P)\£ where 

P 




p 

(c) For some/all A £ p{T^) there exists Dx G £(F, Favi(yl)) such that 

\-A={\~V) 



Id ~Dx 
Id 

Proof. The equivalence of (a) and (b) follows from Engel and Nagel [9j Theorem 
III. 3. 9], while (b) and (c) are equivalent by [Qj Proposition III. 3. 18. (ii)]. Finally, (a) 
and (a') are equivalent since every strongly continuous semigroup is exponentially 
bounded. □ 
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3. Applications 

In this section we will show how our abstract results apply to inhomogeneous 
Cauchy problems as well as to systems with boundary feedback. 

3.1. Inhomogeneous Abstract Cauchy Problems. Consider the inhomoge- 
neous Cauchy problem 



(iACP) 



uit) = Au{t) + f{t) for t > 0, 
u(0) = Uq. 



for a linear operator {A, D{A)) on a Banach space E. For operator splitting methods 
applied to this problem, see Bj0rhus 5 and Ostermann and Schratz [14) . 

A standard method to tackle this problem is to rewrite it as a homogeneous one 
like (|ACPp in the product space £ -.^ E x F(M_|_; E) for the operator matrix 

'^d ^ with diagonal domain D{A) = D{A) x Fi{R+;E). 

Here F{R+;E) denotes a space of i?- valued functions defined on K+ on which the 
left-shift semigroup {L{t))t>o is strongly continuous. Moreover, with domain 
Fi(R^;E) denotes the generator of (i(i))t>o, and (5o(/) ■— /(O) is the point evalu- 
ation at 0. The main choices for F :— F{R.+ ; E) are F — Co(M+; E) which implies 
Fi{m+;E) = C^(R+;E) oi F = L^{R+:E) for which Fi{R+;E) = W^''^{R+;E) 
follows. Then the inhomogeneous equation (jiACPp is equivalent to the abstract 
Cauchy problem 

j ^^u{t) ^ Au{t) fori>0, 
l"(0)-(7). 

For the details we refer to Engel and Nagel 19, Section VI. 7]. Here we only mention 
that in both cases A generates a strongly continuous semigroup on £. In the Co- 
case this easily follows by bounded perturbation (see below) while in the L^-case 
this is shown in W, Proposition VI. 7. 5]. 

3.1.1. Stability in F — Co(R+; E). To see that the operator matrix A is actually a 
generator in case F — Co(IR+; E), note that 

■^0^ ^ with diagonal domain D{Ao) = D{A) x D{-^) 



for D{-^) = Cl{M+;E) generates the semigroup 

Toit) 



(T{t) 
I L{t))^ 



where {T{t))t>o is the semigroup generated by A. Since Sq : F ^ E is bounded, A 
is a bounded perturbation of Aq, hence it is a generator. To get a formula for the 
semigroup (7'(i))t>o generated by A, note that by Proposition 12.11 the semigroup 
(7~(t))t>o rnust be upper triangular, say 



r{t) = 



Ti{t) Ut) 

n{t) 
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By the variation of constants formula (see e.g., Engel and Nagel j9l Section III.l]) 
we obtain 

mi}) = ro(t)(p + 1^ Tit - (^0° ) To(s)(p ds 

If we take / = 0, we get Ti{t) — T{t), and hence for aU / G F we have 

T2{t)f = I T{t- s)f{s) ds. 



Moreover, T3{t) = L{t). 

Now we want to apply the sequential splitting to the problem 



(3.1) 



-f^u{t) = {A, + A2)u{t) + (/i + f2){t) for t > 0, 
u(0) = uo. 



where we have written the inhomogeneity already in a form corresponding to the 
splitting procedure. Namely, choosing a time step h — ^, we first solve the equation 

(^^v(h) = A^vih) + h{h), 
\v{0) = uo, 
then using the result we solve the equation 

if^w{h)^A2wih)+f2{h), 

\wlo) = v{h). 

Setting Uh = w{h), we repeat this procedure n times and call Unh the (sequential) 
split solution corresponding to the equation p.ip . 

Clearly, by the preparations in the beginning of this section, we can reformulate 
p.ip as a homogeneous abstract Cauchy problem 



^w(t) = (^1 + A2) u{t) for t > 0, 

/ Mo \ 

.(0) = (/.) 



on the product space 

£ = Ex F X F, 
for F = Co(K+; E) and the operators 



(3.2) 



A, 


So 


0' 





d 
ds 














''A2 





So 

















d 

ds , 



DiAi) = D{A,) X D{^) X F, 



d N 
ds ' 



for D{-^) — Cq(R+; E). By the above, the semigroups generated by these operators 
take the form 

Tiit) Qi{t) \ ( T2{t) Q Q2{t) 

(3.3) Ti{t) = ( L{t) , T2{t) = / 
1/ \ L{t) 
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where {Ti{t))t>Q and {T2{t))t>o denote the semigroups generated by Ai and A2, 
respectively, {L{t))t>Q is the left-shift on Co(M+;i?), and 

Q,{t)f ^ f Ti{t - s)fr{s) ds for i = 1, 2 and f e F. 
Jo 

Note that with this notation, the sequential sphtting is given by the Trotter product 
formula 

Unh = 7r, {T2{h)Tl{h)T [hj. 

The next result establishes the stability condition (|1.1[) . and hence the convergence 
for the Trotter and Strang product formulas with respect to the splitting A = 
Ai + A2 for Ai and A2 defined by 



Proposition 3.1. Suppose that for some M' > 1 and uj' > one has 



ll(T2(^)Ti(i)r|| <M 

for all t > and n G N. Then there exist M > I and a; G M such that 

(3.4) ||(7^(^)Ti(^))"|| < Afe*" holds for allt>0 and neN. 

Moreover, the product formulas ()1.2|) and ()1.3|) described in Theorem \1.1\ with re- 
spect to the operator splitting A = Ai + A2 for Ai, A2 defined by p.2p hold. 

Proof. Since {Ti{t))t>o and (72(i))t>o are bounded perturbations of diagonal semi- 
groups, the claim follows from Theorem 12.21 and Proposition [221 D 

In a similar way we obtain the following results concerning the weighted splitting. 

Proposition 3.2. Suppose that for some M' > 1 and w' > one has 

||^(Ti(i)r2(i)-Kr2(^)ri(i)r|| <AfV-' 

for all t > and n G N. Then there exist M > I and a; G M such that 

(3-5) \\^{Mim{i) + r2{imi)r\\<Me'- 

holds for all t > and n G N. Moreover, the product formula (|1.5p described in 
Theorem ] 1.1\ with respect to the operator splitting A = Ai + A2 for Ai, A2 defined 
by (1221) holds. 



Proof. The proof follows similarly as the one of Proposition 13.11 from Theorem 
and Proposition 12.61 □ 



Note that the condition w' > is neither a restriction, nor crucial, and was chosen 
only to simplify our calculations in the following subsection. 

3.1.2. Stability in F = L^(R+;i?). Our aim is now to prove that Propositions 13.11 
and 13.21 remain true if we replace the space F = Co(M+;i?) by F = LP(R+;i?) 
for some 1 < p < 00. This is not straightforward since in the Co-case the stabil- 
ity condition p.4p follows by bounded perturbation. However, in the L^'-case the 
perturbation 

-^0 • ^ii^) = W^'f (R+; £;) cF^E 
is unbounded on F and hence it is not guaranteed in general that the off-diagonal 
perturbing term R{t) is 0{t) as i — 0+. Nevertheless, due to a particular additivity 
property of the norm in L^, stability prevails also in this case. To show this, suppose 



PRODUCT FORMULAS FOR OPERATOR MATRICES 



13 



that the conditions of Proposition 13.21 are satisfied. First we group the entries of 
Ti{t), i — 1,2, from p.3p and obtain the 2 x 2-block matrices 





(Tiit) 


Qi{t) 





TiW = 





Lit) 












/ 




(T2{t) 





Q2{t) 


T2it) = 





/ 












Lit) 



Tiit) 




T2it) 




Riit) 
Slit) 



R2it) 

S2it) 



Here, as in the previous case, (L(i))t>o denotes the left-shift semigroup which 
is now defined on the space F — LP(R+;i?) and has generator with domain 
Di-^) = Wi^P(M+;£'). Then from and (HH) we obtain that 



'{T2ih)T,ih))'' (*) (**) 
{T2ih)Tiih)y = Likh) 







Likh)^ 



where 



and 



fe-i 



(*) - E(^2m(/i))^T2(/i)gi(/j)L((fc - J - i)h), 



J=0 



fe-1 



J=0 



In the Co-case Qiih) and Q2ih) were as ft- — ?> 0^ and hence rather crude esti- 
mates for the sums (*) and (**) aheady imphed stabihty. In the present situation 
we have to be more careful and estimate 



fc-i 



^*)|| <5] (r2(ft)ri(ft))^ ■\\Q2ih)L{ik~j-i)h) 



3=0 



k~l 



< M'e"''^''^ ||Q2(/i)i((fc -j - l)h)\ 

3=0 
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Note now that since / G L^, it is also in Ll^^. Using the additivity of the L^-norm 
with respect to the domain of integration we obtain for f (z F 

fc-i 

J2\\Q2{h)L{{k-j-l)h)f\\ 

3=0 

= El /' T2ih - s) [L{ik - J - l)h)f] (s) ds 



3=0 

El 

3=0 



T2{h-s)f{s + {k-j-l)h))ds 



r{3 + l)h rkh 

<M'e-''^Y. ||/(s)||ds = A/'e-'M ||/(s)||ds 

Cj-kh \ p 

J ||/(s)fdsj < Af"e""'=''||/||. 

This imphes that there exists M > 1 and w > so that 

IK**) II < Afe"'^'' and similarly ||(*)|| < Me'^''^ 

hold. 

Remark 3.3. In the same spirit, the stability of the weighted splitting can also be 
established. Since the proof is straightforward and would be only a repetition of 
what we had done so far, we omit it. 

We therefore obtain the following results. 



Proposition 3.4. Propositions \3.1\ and \3.2\ prevail for F — L^(R+: E) and with 
domain D{-^) = W^'P{R+;E). 

Summing up both cases, we have established that the splitting for the inhomo- 
geneous abstract Cauchy problem with a Co(]R+; E) or a LP(R+; E) inhomogeneity 
is stable if the splitting for the associated homogeneous problem is stable. 

3.2. Abstract Boundary Feedback Systems. Let E and dE be Banach spaces 
and let the operators Am : D{Am) C E ^ E, B : D{B) CdE-)-dE,C€ C{E, dE) 
and L : D{Am) — > dE be given. An abstract boundary feedback system is a system 
of two coupled differential equations of the form 

' -l-^u{t)= AraU{t), t>0, 

x{t) = Bx{t) + Cu{t), t > 0, 

Lu{t) = x{t), t > 0, 

^m(O) = uo, a;(0) = xo, 

where the functions u and x are E and i9i?- valued, respectively. We refer to Casarino 
et al. [B] for more details and concrete examples. Now under suitable assumptions 
(see below) such systems can be rewritten as an abstract Cauchy problem (jACPp . 
where the coupling Lu{t) = x{t) of the two equations is coded in the domain of 
the system operator A. To proceed we make as in Casarino et al. [BJ Section 2] the 
following 



(ABFS) 
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Assumption 3.5. (i) A :— Amlka-L generates a semigroup {T{t))t>o on E. 

(ii) L : D{Arn) dE is surjective. 

(iii) (-^Z') : D{A^) ^ E x dE is closed. 

(iv) B generates a semigroup {S{t))t>Q on dE. 

Then by Casarino et al. [6l Lemma 2.2], the following holds. 

Lemma 3.6. If X E piA), then the restriction -/j|kcr(A-A,„) '■ ker(A — Am) dE is 
invertible and its inverse, called Dirichlet operator, 

Dx (^|kcr(A-A„))"' ■■ 9E ^ ker(A - A„,) C E 

is bounded. 

Remark 3.7. We note that condition (iii) in Assumption 13.51 can be replaced by 

(iii') Dx e C{dE, E) exists for all A G p{A) 
which sometimes is easier to veriiy than the closedness of ( £"). 

In order to treat (jABFSp by semigroup methods we define on £ :— E x dE the 
operator matrix 

'Am 0^ 



with domain 



C B 



D{Ac) : = {(D e D{Am) X D{B) -.Lf^x). 

Then by Casarino et al. [6, Section 2] and by Engel and Nagel the system 
(|ABFS[) is well-posed if and only if the operator matrix Ac generates a strongly 
continuous semigroup {Tc{t))t>o on £. Moreover, in this case for every initial value 
("o) ^ D{Ac) the unique solution of (|ABFS|) is given by 

M+3i.->7ri(rc(t)0) GS. 
In order to apply the splitting approach to this problem we first assume that C ~ Q 
and decompose = -^i + A2 for 

0)' D{A^):^{OED{Am)y^dE:Lf^x], 

(3.6) ^ ^ 

M:=(^Q D{A2):=ExD{B). 

Then by Casarino at al. [61 Corollary 2.9] the matrix Ai is the generator of a 
strongly continuous semigroup (7i(i))t>o- Moreover, if A is invertible, then 7i(t) 
is given by 

On the other hand, also A2 is a generator of a strongly continuous semigroup 
('72(t))t>o which can be easily calculated as 

^Id 
S{t)^ 

Now the assumptions (i)-(iii) of Theorem 12.21 and Theorem l2.3l are satisfied if there 
exists K > Q and w G K such that for Ri{t) := {Id — T{t))Do we have that 

(3.7) \\Riit)\\ < Kt ■ e"* for all t > 0. 



T2it) 
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Note that by Proposition 12.61 Lemma [3.81 and [6l Lemma 2.6], condition p.7p is 
equivalent to the assumption 

(3.8) c Favi(^). 

This condition can be characterized by the following result of Desch and Schap- 
pacher [71 Theorem 9]. 

Lemma 3.8. Let the Assumptions [^7S[ fi)-fm) be satisfied. If Dx denotes the 
Dirichlet operator introduced in Lemma [KBl then the following conditions are equiv- 
alent. 

(a) D{A^) C Favi(A). 

(b) ker(A - Am) C Favi(^) for some A e p{A). 

(c) There exist 7 > and Aq £ K such that \\Lx\\ > 7A • for all A > Xq, 
X e ker(A — Am). 

(d) There exist c > and w > such that \\Dx\\ < c ■ A^^ for all X > w. 
Summing up, we obtain the following. 

Corollary 3.9. Let the Assumptions be satisfied. If in addition G p{A), 
C — and D{Arn,) C Fav(A), then the product formulas (jl.2p and (11.31) described 
in Theorem ] 1. 1\ for the decomposition Aq = Ai + A2 and Ai, A2 defined by p.6p 
converge to the semigroup (7o(^))t>o generated by Aq. 

In the next step we add a non-zero feedback operator C G jC{E, dE) to our setting. 
More precisely, we decompose 

(3.9) Ac^Aq+C where ^ - (^c o) ^ 
Then from Proposition 12 .41 we obtain the following result. 

Corollary 3.10. Let the Assumvtions \3.5\ be satisfied and let C € C{dE, E). Then 
the product formulas (II. 2[) . (jl.3p and (11.51) for the Trotter, Strang and weighted split- 
ting with respect to the decomposition (j3.9p converge to the semigroup {Tc{t))t>o 
generated by Ac- 
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